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Abstract 

Langevin equation describing soft modes in the quark-gluon plasma is reformulated on the loop 
space. The Cauchy problem for the resulting loop equation is solved for the case when the nonvan- 
ishing components of the gauge potential correspond to the Cartan generators of the SU (iV)-group 
and are proportional to a constant unit vector in the Cartan subalgebra. The regularized form of 
the loop equation with an arbitrary gauge potential is found, and perturbation theory in powers 
of the 't Hooft coupling is discussed. 
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1 Introduction 



Loop equations in QCD were derived about 20 years ago [|l[] (for a review see ref. [||), and nowadays 
there exist numerous applications of the loop-space techniques to other theories in D > 2. Among 
those, apart from various proposals of string theories obeying loop equations |3|, §|, it is worth 
mentioning applications to supersymmetric gauge theories p, [J| , gravity M , and turbulence |7]. Q . 
In the present letter, we shall derive and analyse a novel loop equation describing the dynamics 
of soft gauge fields in quark-gluon plasma. It is based on the following Langevin-type equation H 

1 d t A i = -D j F ji + g( i , (1) 
where the bilocal correlator of Gaussian noises reads 

(ff(s)C j(y)) = 2T 1 5 ij 6"»6(x - y). (2) 
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In these equations, i,j = 1, 2, 3, a = 1, . . . ,N 2 — 1, 7 is the so-called color conductivity, whose 



value in the pure SU (iV)-theory under study is |10 



4vrT 

T is the temperature, A{ = igAft a is an anti-Hermitean matrix with tr t a t b = 5 ab , Di = <9j + A4, 
and Fij = diAj — djAi + [A iy Af\ (here we have adapted the notations of refs. (T|, 

Let us briefly quote the qualitative derivation of eqs. (Jj])-(@) presented in ref. [[□]]. It is based on 
the idea that soft modes are effectively classical and therefore can be described by the non-Abelian 
analogue of the Maxwell equation D x B = D t E + Jhard- Here, B = D x A and all the covariant 
derivatives involve only soft degrees of freedom, while the color current Jhard describes hard degrees 
of freedom. Owing to the fact that plasmas are conductors, this current can be expressed in terms 
of E as Jhard = lE. Assuming that the characteristic temporal scale of soft modes is much larger 
than their spatial scale (that will be discussed later) and adapting the gauge A = 0, we arrive 
at the following equation: ^d t A = —D x B. This equation is dissipative, and the dissipation 
is a consequence of interaction of the soft modes with the hard ones. However, hard modes not 
only acquire the energy from the soft modes, but also serve as a source of thermal noise (. The 
condition necessary to reach the equilibrium distribution of soft modes is the balance between 
their excitation by the thermal noise and dissipative decay. For a general Langevin equation of 
the type ^d t q = — V $V(q ) + one can show JE2] that the correct equilibrium distribution e~ v ^ T 
is reproduced by the Gaussian noise, whose correlation functions read 



/2fc+l \ / 2k 



m)Q(t')} = 2^5^ -if), ( n cm ) = o, ( n cm > = e n <& &) & &)> • 

\ a=l I \a=l I possible pair pairs 

combinations 

These are the qualitative considerations that lead to the effective theory for soft modes, described 
by eqs. ©-©. 

Let us also discuss the characteristic spatial and temporal scales of soft fields. Since we are 
interested in soft (or nonperturbative) fields, both terms in the covariant derivative associated 
with such fields, should be equally important. This means the inequality A > 0(1/ gR), where 
A is the amplitude of a typical nonperturbative fluctuation of the spatial size R. The energy of 
such a fluctuation is thus E ~ RA 2 > 0(1/ g 2 R). For the probability of fluctuation not to be 

suppressed as e~ E ^ T ~ e 3 2rt , one needs R > 0(1/ g 2 T). Therefore, the characteristic spatial 
scale of soft fields is 

R ~ 4p (4) 
g 2 T 

Next, from eq. ([I]), we have the following estimate involving the temporal scale r of soft fields: 
7r _1 y4 ~ R~ 2 A, so that r ~ 7.R 2 ~ 7/ (g 4 T 2 ). By virtue of eq. (J3|), we finally have 



In particular, this yields the estimates 



-A fl„ „-lV //I T . . „2> -1 



7 a V a 1 ' R 



g*ttng- 1 ) « 1, — - ~ g A In g~ l « 1, 
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where the inequalities are obviously valid for small enough g. These inequalities justify the omis- 
sion of the term D t E in the equation D x B = D t E + jE above. 

The letter is organized as follows. In the next Section, we shall present the derivation of the 
large- N loop equation in quark-gluon plasma. In Section 3, the Cauchy problem for this equation 
will be solved at a certain class of fields. In Section 4, the regularized version of the equation will 
be derived, and the perturbative expansion will be briefly discussed. Finally, the main results of 
the letter will be presented in Summary. 



2 Derivation of the loop equation in quark-gluon plasma 

Similarly to the theory of turbulence [0], the Wilson loop, 



W = W[C,t] = ^ (tiV exp dxiA^t)^ 



is a functional of the 3D contour C and a function of time. It it worth emphasizing that the 
average here is implied in the sense of the average over the Gaussian noise, 



(O) 



rather than over the Yang-Mills action. 

Differentiating W with respect to time and using eq. ([!]) we get 



JVQOexp 


oo 

-^-Jd'xJdtQ 2 

' 


J VQ exp 


oo 

-^Jd'xfdtQ 2 

' 



(6) 



(7$ + A) W = jL j dx t tr (C«(x, t)U xx ) . 
c 

(J)6xi((r) * s tne functional Laplacian [] and 
Uxy = V exp J dziAi(z,t) 

\Pxy / 

is the parallel transporter factor along the piece C xy of the contour C which starts at the point y 
and ends up at the point x. The correlator tr (C»(x, t)?7xx) on the r.h.s. of eq. (|>D can further be 
rewritten as follows |: 



(7) 



Taking into account that \V1 



1 Rigorously speaking, this form of the functional Laplacian is admissible only when it acts onto functionals like 
W , which do not have marked points. 

2 Although we do not expect any confusion between the space and matrix indices, i,j, note that the former are 
always denoted by the lower case, while the latter are denoted by the upper one. 
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<5<?(x,f) 2 7 " 

we arrive at the expression gT(t a ) %: > (^ SA J^ ^y ^ X x)- ^ can evaluated by virtue of the formulae 



and 



(t a )^'(t a ) fc ' = 5 a <^ - — 5* J '(5 W . (9) 

Using the large-iV factorization of the resulting two-loop average, we eventually arrive at the 
following large- N loop equation equivalent to the original eq. ([[]): 

(7$ + A) W = XT <j) dxi <j> d yi 5(x - y)W[C xy , t}W[C yx , t\. (10) 

c* c 

Here, A = g 2 N is the 't Hooft coupling constant which remains finite in the large- N limit, and the 
contours C xy , C yx are always closed due to the presence of the 5-function. Note that functional 
form of the r.h.s. of eq. (fLOD is the same as that of the usual loop equation jjj, while the l.h.s. is 
different due to the presence of the temporal derivative. 



3 Solution of the Cauchy problem for a certain class of 
fields 

The Cauchy problem for eq. ( |T0D with the initial condition W[C, t = 0] = Wo[C] can be solved 
for the case when the only nonvanishing color components of the gauge potential correspond to 
the Cartan generators H = (H 1 , ■ ■ ■ , H N ~ 1 ^j of the group SU(N), i.e., A4 = AjH, and moreover 
Aj = mpi. Here, n is a constant unit vector in the Cartan subalgebra, and ifti(xi, t) is an arbitrary 
3D vector-function. Introducing the notations <£> xy = f c dziipi(z,t), (^yx = /c dziipi(z, t) and 
making use of the obvious formula e lt/mH = cos ip+inH sin (p (which stems from the ortho normality 
of roots |TI|, HiHj = Sij), we have in the large- N limit under study 

W[C xy ,t]W[C yx ,t] = (e iVxy ) (e^ yx ) ~ ^ e i(Vxy+ ^ x) ) = W. 

Denoting further V = InW and using the fact the the functional Laplacian obeys the Leibnitz 
rule, we can rewrite eq. ([TOD as 

frd t + A)V[x,t] = J[x], (11) 

where J[x] = XT § c dxi § c dyiS{x — y) and from now on we shall use an alternative notation 
0[C] = O [x], so that x = x(a), < a < 1 is an element of the loop space. 

In order to invert the operator standing on the l.h.s. of eq. (|TTJ) (that would lead to the solution 



of the Cauchy problem for this equation) we shall make use of the method proposed in ref. [14 



Note that this method has been applied in ref. M in order to solve the Cauchy problem for the 
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loop equation in turbulence, and here we shall use the same strategy. In ref. [13], it has been 
shown that for an arbitrary functional O [x] defined on the loop space, the following equation 
holds 



\X + 

In this equation, the smeared Laplacian A^ has the form 



d , , 



X 



'12) 



A< G > = J da v.p. j da'G (a - a') 



5x, (a') Sx, (a] 



+ A, 



(13) 



o o 

where v.p. / da' = Jq CT ~° da' + /J +0 da', and the average over loops is defined as follows 



O 



(14) 



e(o)=l*(i) 

Here, G (cr — cr') is a certain smearing function and 

5 = - / rfa | dof (a) CT 1 (a - a') |* (<t') 
o o 

with G _1 standing for the inverse operator. Note that the first term on the r.h.s. of eq. (fl3D is 
an operator of the second order (i.e., it does not obey the Leibnitz rule) and is reparametrization 
noninvariant. 

In particular, for 



(15) 



G (a - a') 



£<1, (16) 
. For G (a — a') defined by eq. (|16|), 



the action (|i~5f) becomes local: S = \ J da e (d (T ^') + -£ 

o L v ' 6 

the contribution of the first term on the r.h.s. of eq. ( |T3"D is of order e for smooth contours, 
and therefore when e — ^ it vanishes, A( G ) tends to the standard functional Laplacian A, and 
reparametrization invariance becomes restored. In general, the smearing function G(a) obeys the 
relations G(a ^ 0) — > at e — > 0, (7(0) = 1 and can be obtained from the one given by eq. (|i"6|) 
upon a reparametrization. Equation fll2"D is then a consequence of the equation of motion 



da'G (a - a') 



'SO 



(17) 





which follows from eqs. (|14] ) and (|15"1). 

Performing a shift of the contour C, a? — > x + \/A£ in eq. (|TT]) and using eq. fll^f), we arrive at 
the following equation: 



' <9 N 



r 



X 



j 



X 



The solution to the Cauchy problem for this equation obviously reads 



5 



V 



x 



V 



x + JA--Z,0 



^ + 7 
£ 7 o 



2 + JA--(t-T)£ 



Taking further into account that for nonperturbative fields under study, the minimal area asso- 
ciated with the contour C, S m i n , obeys the inequality S m i n ^> R 2 and recalling eqs. we 
obtain the following estimate: 



A(V 



x 



T 



T 
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7 In g 1 An 



0.72. 



Using this smallness and eq. flT7|), it is straightforward to get the following formula: 

2N 



V 



x + a A 



2t 

7 



1 _ I A (G) + O 

7 




after which 



Obviously, an analogous expansion can be derived for ^ J x + y'A — ^(t — r)£ 

one can send the parameter e to zero and take into account that A J [x] = 0. Setting finally the 
parameter A to be equal zero (that removes the ^-dependence completely) and passing from V 
back to W, we arrive at the following solution to the Cauchy problem for eq. (|TT|): 



W[x,t] = W [xjexp 
In particular, by virtue of the formula [2] 



7 V' 1 ' 



awq g; 







(18) 



A <j) dxi <j> dyiD (x - y) = -2 <j> dxi I dyiS(x - y) 

c c c a 



(19) 



where -Dq( x ) = (47r|x|) is the massless propagator, it is straightforward to see that at Wq [x\ 



XT 



exp ^— f c dxi § c dyiD (x. — y)j, W[x, t] = Wo[x]exp O ^( 7 / . j j ■ It is also worth noting 

once more that the obtained solution flT8f) is valid only at t < 0(t), rather than at an arbitrarily 
large t. 



4 Regularized loop equation and perturbative analysis 



In order to regularize eq. (|T0"|), one should use the regularized version of eq. ([[]). It is analogous 
to the one proposed for the method of stochastic quantization in ref. |15| and reads 



7 ^(x,t) = -D*F%(x,t)+g j d*yR(DX y C(y,t). 



ab 



Here, R{D 2 )'^ y = ^e^J 5(x — y) is the regularizing function with A standing for the UV cutoff. 
An analogous method of regularization of the usual loop equation has been used in ref. [[R| , where 



3 This fact makes the expansion for the current J [x] exact at the zeroth order of the parameter 
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the equality of —A to the field-space Laplacian (implied in the weak sense) has been employed. 
Here, we shall perform the regularization directly, without the use of correspondence between the 
two Laplacians. One has 

A- 2 

-if drv 2 (r) 

R(D 2 ) a ^ y = J Vre ° tr 

r(A-2)=y 
r(0)=x 

where the parallel transportes are defined as in eq. (^), but along the regulator paths r xy and r yx , 
whose typical length is of the order A -1 . By virtue of this representation, we obtain the following 
analogue of eq. (|6|): 

(jd t + A) W = 



t a U(r xy )t b U(r yx ) , 



A" 2 

n r r r -5 f dTi- 2 {r) 

= jjfdxij d 3 y J Vre ffiC)"^ (C b (y, t)U lm (r xy )U nk (r yx )U^ x ) . 

C r(A"2)=y 
r(0)=x 

Having in mind the eventual limit A — > 00, we may use eq. (^) for the evaluation of the latter 
average. Indeed, in this case, the use of the regularized expression instead of the 5-function in 
eq. (H) would lead to the excess of accuracy in the (A — > oo)-limit (i.e, to the leading order in 
1/A under study, the resulting nonlocal expression will anyway be reduced to the local one). By 
virtue of eq. (Bl), we then obtain in the large- N limit: 



A" 



( 7 ft + A) W = 9 — f dxi J d 3 y J Vre / j^^U^ (r xy )U n \r yx )U^ \ . 

C r(A~ 2 )=y \ ' • 

r(0)=x 

In the evaluation of the last variational derivative, we may take into account that the paths r xy 
and r yx are infinitesimal, i.e., to the leading order in 1/A it is enough to differentiate only £/^ x . 
Integrating then finally over d 3 y by virtue of the resulting 5-function, we arrive at the following 
regularized equation: 

A" 2 

r r r -a / dTt ' 2 ^ 

(jd t + A)W = XTfdx i fdy i J Vre W[C xy r yx ,t]W[C yx r xy ,t]. (20) 

C C r(A"2)=y 

r(0)=x 

Clearly, in the limit A — ► 00, it recovers eq. (p!Q|) in the same way, as the regularized equation 
derived in ref. J16 recovers the usual large-A^ QCD loop equation. Owing to the closeness of the 
contours C xy r yx and C yx r xy , eq. (|20D is formulated entirely on the loop space. 

Let us now turn ourselves to the perturbative analysis of the formal solution to the Cauchy 
problem for eq. (|I0| ). This solution obviously reads 



7 



1 



W [x, t] = - I dre^r A J [x, r] + e~$ A W [x] , (21) 
7 J 
' o 

where J[x,t] is defined by the r.h.s. of eq. (|T0D. Let us assume that d t W vanishes at t — > 
and justify this assumption later perturbatively. Then in the limit t — > 0, eq. (|ToD takes the 
form of the usual 3D loop equation (with the only difference that the dimensionful 3D 't Hooft 
coupling is replaced by AT), which is obeyed by W . Owing to this fact, to the order A 1 , W can 
be represented as [cf. e.g. ref. M and eq. (O)] f\ 

XT r r 

W a = l- — j dxij d yi D (x -y) + (A 2 ) . (22) 

c c 

Substituting further this expansion together with the Ansatz W = + AW^ 1 ) into eq. ([H]), we 
immediately find = 1, while for we have 

W {1) = — j> dxij dyi5(x-y)- — M - -Aj j> dx { j> %D (x-y) = —^j d%i / %-°o(x-y). 

^ c c ^ ^ ' c c c c 

In the derivation of this expression, apart from the formula A § c dxi § c dyi5(x — y) =0, mentioned 



in the previous Section, we have again used eq. (191). Thus, we see that W = W up to the first 



order of the expansion in A. In particular, this fact justifies our initial assumption lim^Q dtW = 
at least up to this order of A. The t-dependence of W appears in the next orders of perturbation 
theory. This is because the action of A™, (n = 2, . . . , n max , where n max is determined by the order 
of perturbation theory) onto the respective terms of the expansion of Wq gives then a nonvanishing 
result. 



5 Summary 

In the present letter, we have derived and analysed the large- N loop equation in quark- gluon 
plasma. This equation is based on the Langevin equation which describes the dynamics of soft 
degrees of freedom. The regularized form of the loop equation, formulated entirely on the loop 
space, has also been obtained. It is derivable directly from the regularized version of the original 
Langevin equation. By making use of the method of smearing of the loop-space Laplacian, we 
have further solved the Cauchy problem for the unregularized equation in some particular case. 
Namely, this has been done for the gauge potential of the form Ai = nH^ i; where H's are the 
Cartan generators of the SU(N)-group, n is a (N — 1) -dimensional constant unit vector, and 
ijji(~x,t) is an arbitrary 3D vector-function. The obtained solution is approximate and valid up to 

terms of the order O ) ^ , which are small for the characteristic temporal and spatial scales 

of soft fields. We have also briefly addressed the perturbative expansion of the formal solution 
to the Cauchy problem in the general case (i.e., for an arbitrary gauge potential). Clearly, to 
perform this analysis, one first needs to fix the perturbative expansion for the initial value of 

4 Clearly, another assumptions imposed on the limit lim, n d t W correspond to the choices of Wq which yield 
the perturbative expansions different from the one of eq. (|22| ) ■ 



s 



the Wilson loop. In order to extract such an expansion from the obtained loop equation, it is 
necessary to postulate a certain spatial dependence for lim^Q d t W . In the simplest case, when this 
limit is set to zero, the loop equation for the initial value of the Wilson loop takes the functional 
form of the usual large-A QCD loop equation. The latter is known to be satisfied by the 
perturbative expansion of the Wilson loop, computed in large-A QCD. Using this expansion, it 
has been checked that for the case lim t _^ d t W = under study and to the first order in the 
't Hooft coupling, the solution to the Cauchy problem is equal to the initial value of the Wilson 
loop at any times. (In particular, this fact justifies the postulate lim^ d t W = itself, within this 
order of perturbation theory.) However, even in such a case, the time-dependence of the solution 
does appear in the higher orders of perturbation theory. 
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